Minimal measurements of the gate fidelity of a qudit map 
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We obtain a simple formula for the average gate fidelity of a linear map acting on qudits. It is 
given in terms of minimal sets of pure state preparations alone, which may be interesting from the 
experimental point of view. These preparations can be seen as the outcomes of certain minimal pos- 
itive operator valued measures. The connection of our results with these generalized measurements 
is briefly discussed. 
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The interaction with the environment has a significant 
negative impact on any physical implementation of a uni- 
tary gate or a quantum channel. Dccoherence, turning 
pure states into mixed states, creeps up and unitarity is 
lost in the Hilbert subspace of the signal states. The char- 
acterization of the quality of such implementations is, 
hence, of utmost importance for quantum computation 
from both experimental and theoretical point of view |Q . 

A physical quantum gate or channel is best described 
by a trace-preserving linear map or superoperator £, 
which is assumed to be an approximation of a unitary 
operator U characterizing the quantum gate (U = 1 for 
a quantum channel). The average gate fidelity 

F(£, U) = J dip tr [Up^£{p^)] , (1) 

where dip is the invariant Haar measure on the space 
of pure states \ip) (i.e., dtp = dtp' if \ip') = U\ip)) and 
Pf = I^X^Ij is recognized as a convenient figure of merit 
and is widely used to assess the quality of a quantum 
gate or channel. 

For gates acting on qubits, Bowdrey et al. j|] have re- 
cently derived a simple formula for F(£,U), from which 
they obtain a convenient way of actually measuring the 
gate fidelity in a laboratory. It amounts to replacing 
the integral in ([!]) by a sum over a finite number of 
pure states whose Bloch vectors point at the vertices of 
an octahedron or a tetrahedron inscribed in the Bloch 
sphere. Nielsen || has obtained a similar formula for qu- 
dits (i. e., quantum states belonging to a d-dimensional 
Hilbert space) in terms of a finite set of unitary opera- 
tions Uj orthogonal with respect to the Hilbert-Schmidt 
inner product, namely, such that tr (UjUj) = dSij. His 
derivation is based on Horodeckis' formula connecting 
F(£,l) = F(£) with the entanglement fidelity 0, for 
which he also provides a simple proof. Our aim is to 
give a more natural generalization of Bowdrey's formula 
so that only measurements over a (minimum) number of 
pure states need to be performed to verify F experimen- 
tally. Interestingly enough, the solution of this problem is 
not far removed from that of obtaining minimal positive 
operator valued measurements (POVMs) in the context 



of optimal communication of directions through a quan- 
tum (spin) channel, which has received much attention 
over the last few years || |[ f?J. 

In the first part of this letter we derive a general ex- 
pression of F in terms of the SU(d) group generators. In 
the second part we write F as an average of measure- 
ments over a finite and minimal number of pure state 
preparations, which may be experimentally relevant. 

We first notice that due to the invariance of dip one has 
F(£,U) = F(£',t) = F(£'), where £'{p) = £{WpU). 
Hence we only need to consider the simpler form F{£ ) 
without any loss of generality. Furthermore, the uniform 
measure dip can be effectively realized in the following 
way: J dipT{ip) = J dll ^(Uipo), where T is any func- 
tion of \ip), \ipo) is a fixed reference state, and dll is the 
Haar measure of SU(d) normalized such that J dll = 1. 
This ensures that J dipT(ip) — J dipT(U'ip) for any 
U' € SU(d). Obviously, not all the d 2 — 1 parameters 
involved in dU are physically significant. E.g., for qubits 
(d = 2), dip = dn is the uniform measure on the 2-sphere 
S 2 , which can be parametrized by the Euler angles a and 
(3. If \iJjq) is an eigenstate of the Pauli matrix a z , any 
function of \ip) is independent of the third Euler angle 7. 
Hence, J dip T = [J d-y/{2ir) = 1] X / dnT = J dXJ T . 
For qudits, dn must be replaced by the invariant mea- 
sure of SU(d)/[SU(rf — 1) x U(l)], since any reference 
state 1-00) is now invariant under SU(d— 1) x U(l). Note 
that the number of parameters match, as qudits depend 
on 2(d — 1) real variables. With all this we finally have 

F(£) = JdUtY [Up Q U^£(Up U^] , (2) 

where we have defined po = |'0o)('0o|- 

There exists a (gP — l)-dimensional unitary vector with 
components rig, a = 1, 2, ... d 2 — 1, such that the density 
matrix po can be written as 

Po = -; + k d n^T a = -1 + k d n ■ f, (3) 

where kd = y/2(d — l)/d and {T a } are the (hermitian 
and traceless) generators of SU(d), normalized so that 
tr (T a T b ) = S ab /2. They can be chosen as T a = X a /2, 
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where A a are a generalization of the Gell-Mann matrices 
of SU(3). Throughout this letter a sum over repeated 
latin indexes is understood. It is straightforward to ob- 
tain 



n = — tr(poT'). 
kd 

We now recall the well known relation 
UT a tf = {AdU) a b T b , 



(4) 



(5) 



where Ad stands for the adjoint representation of SU(rf). 
We also recall the orthogonality of the irreducible repre- 
sentations of compact groups, which in the present case 
implies 



dU(AdU) a b = 0, 
dU{AdU) a b {AdU) c d = 



6acS M 



d 2 -l 

Using that £ is linear and trace-preserving one gets 

d 2 -l 



F{£) 



d(d+l) 



- tr [Ta£(T a ) 



This is the generalization of Bowdrey's formula 



F{£) 



1 1 X 

9 + ^ E 



tr 



2 t \2 



(6) 
(7) 

(8) 
(9) 



for qubits. That concludes the first part of this letter. 

The key ingredient of the derivation above is the or- 
thogonality relation of the group representations, exem- 
plified by Eqs. |6| and [7| We will show that it is possible 
to find a discrete version of these equations. Namely, one 
can find a finite set of SU(d) elements {U r } and positive 
constants {c r } such that 



J2cr(AdU r )a b = 



Y,Cr(AdU r ) a b (AdU r ) c d = 



Sac* 



bd 



d 2 -l 



(10) 



(11) 



With this we can reverse the steps going from (||) to (||) 
using the relations ( flO|) and ( |TT| ) instead of their contin- 
uous version (0) and m), and obtain 



F{£) = ^c r tr [p r £(p r )] 



(12) 



This equation has a very convenient form which allows 
setting up experimental tests to determine the fidelity of 
a gate or channel, as will be discussed below. 

Let us briefly discuss the solutions to Eqs. [h] and [fl]. 
The idea is to generalize the concept of a finite set of 
isotropically distributed unit vectors introduced in B and 



adapt it to the problem at hand. A sufficient condi- 
tion for ( |l2| ) to hold can be obtained by contracting ( |l0| ) 
and (O) with n^. . If we define 



n%(AdU r ) a b = n b r , 
conditions (Eof) and (O) lead to 



= 0; 



d 2 - 1' 



(13) 



(14) 



In this sense, we may qualify the set {n r } as isotropically 
distributed (as far as the adjoint representation SU(d) is 
concerned). We have traded the problem of finding {U r } 
for that of finding {n r }. Note, however that the set of ma- 
trices {AdU} is a proper subgroup of SO(c? 2 — 1). Hence, 
not any vector on the sphere S rf ~ 2 is admissible. We 
will come back to this issue below. The relations dl(j), 
( |TT| ) and ( pT| ) also appear in the rather different context: 
the construction of finite positive operator valued mea- 
surements which are optimal for communicating a direc- 
tion |(| ^ H . Leaving aside an overall trivial normaliza- 
tion (from Eq. 14 it follows that ^2 r c r = 1 whereas in || 
J2 r Cr = d), the results in those papers can be readily 
used here. In particular, it is proved in || that solutions 
of (|l4|) exist and the minimal one is given by a set of d 2 
vectors pointing at the vertices of a regular hypertetra- 
hedron or, more properly, (d 2 — l)-simplex, inscribed on 
gd -2 Tj^jg hypertetrahedron is defined by the condition 



d 2 - r 



r ^ s, 



(15) 



and the exact overall orientation has to be chosen so that 
all vectors ri r are of the form (|l3|). For this hypertetra- 
hedron all the coefficients c r are equal: in our notation 
Cr = 1/d 2 , r = 1,2,...,g? 2 . A explicit form of n r for 
SU(3) can be found in S, where also the general case 
is briefly discussed. The solution is more conveniently 
expressed in terms of states |?/v) such that 



\l/j r )(lp r \ = p r = UrPolll = - + k d ri r ■ f. 



Then ( |T5| ) translates into 

|^r|^)| 2 = 



d+1 



r 7^ s. 



(16) 



(17) 



Since this equation is a condition on states, Eq [13| is au- 
tomatically satisfied by the corresponding Bloch vectors. 
For SU(3) a solution of dT7| ) has the simple form 



KM - T2 



1 

3 2(r-l)7ri/3 




r = 1,2,3. 



(18) 



The remaining six states are obtained by applying cyclic 
permutations to the components of \1p1-3} (a different 
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choice of states is given in || ) . From a experimental point 
of view, the states ([TJ]) have a very appealing form; each 
of them involves a linear combination of only two states 
of the computational basis. If the qutrit is implemented 
by say three atomic levels, only two levels need to be 
manipulated to prepare each one of the p r . 

In the SU(2) case, a solution in terms of n r is more 
simple and transparent, mainly because SU(2) is isomor- 
phic to SO (3), which makes ( |l3| ) trivially satisfied by any 
vector of S 2 . A compact solution is given by 



ni 



>.u>, 



n 2 



1,-1,1), 



(19) 



and ris, 7T4 are again obtained by applying cyclic permu- 
tations to the components of rii. 

For the hypertetrahedra discussed above, Eq. [l2| can 
be cast as 



F{£) 



1 d 

]2^2 tr \Pr £ iPr)] 



d 2 



(20) 



r=l 



which is the SU(d) generalization of Bowdrey's SU(2) for- 
mula. This equation is our main result and provides a 
remarkably simple procedure for measuring F(£). One 
just has to average the fidelities for d 2 isotropically dis- 
tributed pure states. Notice that all state preparations 
p r have the same weight in the average, thus reducing 
systematic errors. In some sense, Eq. |20| can be regarded 
as the average survival rate of the states {p r } in a quan- 
tum channel characterized by the linear map £. There 
is a alternative way of writing (pC|) which may provide 
further insight. It is straightforward to verify that the 
projectors O r = p r jd are the complete set of positive 
operators (i.e., ^20 r — 1) of a minimal POVM. Thus, 
the states p r are just the preparations produced by the 
device characterized by {Or). We can write (poh as 



_ 1 d * 



We readily see that the same device {O r } could be used 
for the preparations of the pure states {pr} as well as for 
the measurements over {£(p r )}. 



Although we have presented the results for the mini- 
mal sets of p r , for practical reasons, one may wish to use 
a larger number of states. This possibility is easily imple- 
mented in this framework, as conditions ( |l3|) and (14) are 
not specific of minimal sets but entirely general. E.g., for 
qubits one can find a set of six states whose Bloch vectors 
point at the vertices of a regular octahedron || ||, Qj). In 
this case there is a simple setting for preparing the states, 
just three Stern-Gerlach's oriented along the three or- 
thogonal directions. 
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